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Abstract 

A Carleman function is a special fundamental solution with a large parameter 
for the Laplace operator and gives a formula to calculate the value of the solution 
of the Cauchy problem in a domain for the Laplace equation. The probe method 
applied to an inverse boundary value problem for the Laplace equation in a bounded 
domain is based on the existence of a special sequence of harmonic functions which 
is called a needle sequence. The needle sequence blows up on a special curve which 
connects a given point inside the domain with a point on the boundary of the domain 
and is convergent locally outside the curve. The sequence yields a reconstruction 
formula of unknown discontinuity, such as cavity, inclusion in a given medium from 
the Dirichlet-to-Neumann map. In this paper an explicit needle sequence in three 
dimensions is given in a closed-form. It is an application of a Carleman func- 
tion introduced by Yarmukhamedov. Furthermore, an explicit needle sequence in 
the probe method applied to the reduction of inverse obstacle scattering problems 
with an arbitrary fixed wave number to inverse boundary value problems for the 
Helmholtz equation is also given. 
AMS: 35R30, 35R25, 35J05, 33E12, 35C05 
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1 Introduction 

In [24] Isakov considered the uniqueness issue of the problem for recovering the unknown 
inclusions in a given electric conductive medium from infinitely many boundary data. This 
is a special, however, very important version of the Calderon problem [5] and is formulated 
as an inverse boundary value problem for an elliptic partial differential equation in which 
the boundary data are given by the associated Dirichlet-to-Neumann map. Therein he 
established a uniqueness theorem. The points of his paper are: a contradiction argument; 
orthogonality identities deduced by denying the conclusion; singular solutions; unique 
continuation theorem; Runge approximation property; uniqueness of the solution of the 
Cauchy problem for elliptic equations. See [12] for these points. 
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After his work for almost ten years nobody mentioned the reconstruction formula of 
unknown inclusions. However, finding such a formula is a natural attempt in studying 
inverse problems. So the author was wondering why people after him did not immediately 
start their study with seeking such a possibility. However, the fact is: after ten years in 
[10] the author found the full reconstruction formula of unknown inclusions in a conductive 
medium with a constant conductivity. The author named the method developed therein as 
the probe method. It is based on the existence of a special sequence of harmonic functions 
which is called a needle sequence and plays a role of a probe needle. The existence has 
been ensured by using the Runge approximation property which can be proved by showing 
the denseness of the ranges of infinitely many linear integral operators. Therefore the 
construction of the needle sequence is not explicit. 

Succeeding to the probe method, the author introduced another method which he calls 
the enclosure method [13]. This method gave us a different way of using the exponential 
solutions for elliptic equations from Calderon 's way [5] and yielded an explicit extraction 
formula of the convex hull of the unknown inclusions. This opened Pandora's box since if 
one replaces the exponential solutions with any other solutions having a similar property, 
then one can get several information about unknown inclusions. The paper [17] is just a 
result in which, solutions coming from Mittag-Leffler's function 

oo 

instead of the exponential solutions were employed. The result yielded more than the 
convex hull and the numerical implementation has been done in [23] . However, this result 
is restricted to the two dimensional case. 

In this paper we return to the probe method. Our aim is to make this method 
completely explicit. The idea came from reconsidering the role of a Carleman function 
which is a special fundamental solution with a large parameter for the Laplace operator 
and gives a formula to calculate the value of the solution of the Cauchy problem in a 
domain for the Laplace equation. Roughly speaking, we say that a function $(y,x,r) 
depending on a large parameter r > is called a Carleman function for a domain Q 
and the portion T of dQ if, it satisfies the equation Am + S(y — x) = in Q; for each 
fixed x e n <E>( • , x, r) and d/dv§( • , x, r) on dQ \ F vanish as r — > oo. This function 
gives a representation formula of any solution of the Cauchy problem for the Laplace 
equation in Q by using only the Cauchy data on T and yields a natural regularization 
for the numerical computation of the solution from a noisy inaccurate Cauchy data on 
T. The Cauchy problem for the Laplace equation is a fundamental and important ill- 
posed problem appearing in mathematical sciences, engineering and medicine. Therefore 
it is quite important to seek an explicit Carleman function in several domains. For 
this problem Yarmukhmedov [32] gave a very interesting Carleman function for a special 
domain in three dimensions which is a special version of his fundamental solutions [30] 
for the Laplace operator. 

In this paper, using his Carleman function, we give an explicit needle sequence in three 
dimensions. This means that it is written in a closed-form by using an integral involving 
a special function (Mittag-Leffler's function). We never solve integral equations in any 
sense. 

We can summarize the conclusion of this paper in one sentence. 
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• the regular part of a special Carleman function which has been introduced by 
Yarmukhamedov yields the desired needle sequence. 

This makes an explicit link between the probe method and a Carleman function. 

A brief outline of this paper is as follows. In Section 2 we explain what is the probe 
method and the role of the needle sequence. For the purpose a simple inverse boundary 
value problem for the Laplace equation is considered. Needless to say, a reader who is 
familiar with the probe method can skip this section. The main idea is described in Section 
3 by considering a two dimensional case. In Section 4 we introduce Yarumukhamedov's 
fundamental solution for the Laplace operator in three dimensions and show that the 
regular part of a special version of his fundamental solution yields an explicit needle 
sequence. In Section 5 we consider how to construct the needle sequence for the Helmholtz 
equation. We show that the Vekua transform of the regular part of the special fundamental 
solution for the Laplace operator yields the desired needle sequence. In Appendix, for 
reader's convenience, we give a direct proof of a theorem established by Yarmukhamedov. 

2 Probe method and needle sequence 

This section is devoted to a reader who is not familiar with the probe method. Recently 
in [18, 20, 21] the author reformulated the probe method and further investigated the 
method itself. It became much simpler than the previous formulation in [10]. In this 
section, for simplicity of description we restrict ourself to the case when the unknown 
discontinuity is coming from cavities. 

Let f2 be a bounded domain of R m (m = 2, 3) with a Lipschitz boundary. Let D be a 
bounded open set of Q such that: D is contained in f2; dD is Lipschitz; Q\D is connected. 
Given / e ^^(dVL) let u E H 1 (il\D) be the unique weak solution of the problem 



The map Ad : / i — > du/dv\Qn is called the Dirichlet-to-Neumann map. We set A# = A 
in the case when D is empty. 

In short, the probe method is a method of probing inside given material by monitoring 
the behaviour of the sequence of the energy gap 



for a specially chosen sequence {v n } of solutions of the governing equation for the back- 
ground medium (D = 0) which play a role of probe needle. 

The method starts with introducing 
Definition 2.1. Given a point x G Vt we say that a non self-intersecting piecewise linear 
curve a in Q is a needle with tip at x if a connects a point on dQ with x and other points 
of a are contained in Q. 

We denote by N x the set of all needles with tip at x. 



Au = in ft \ D, 




u = f on dQ. 
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Let b be a nonzero vector in R m . Given x G R m , p > and 9 g]0, 7r[ set 
0/2) = G R m | (y - x) ■ b > \y - x\ \b\ cos(#/2)} 

and 

B p (x) = {ye R m \\y-x\ < p}. 

A set having the form 

V = B p (x)nC x (b,9/2) 
for some p, b, 9 and x is called a finite cone with vertex at x. 

Let G(y) be a solution of the Laplace equation in R m \ {0} such that, for any finite cone 
V with vertex at 

/ \VG{y)\ 2 dy = oc. (2.1) 

Hereafter we fix this G. 

For the new formulation of the probe method we need the following. 
Definition 2.2. Let a G N x . We call the sequence £ = {v n } of iJ 1 (f2) solutions of the 
Laplace equation a needle sequence for (x, a) if it satisfies, for any compact set K of R m 
with K C Q \ a 

\™j\v n { -)-G(.-x) \\ LHK) + || VK( • ) - G( ■ - x)}\\ LHK) ) = 0. 

The existence of the needle sequence is a consequence of the Runge approximation prop- 
erty for the Laplace equation. 

In [18] we clarified the behaviour of the needle sequence on the needle as n — > oo. 
The two lemmas given below are the core of the new formulation of the probe method 

Lemma 2.1. Let x £ Q be an arbitrary point and a be a needle with tip at x. Let £ = {v n } 
be an arbitrary needle sequence for (x,a). Then, for any finite cone V with vertex at x 
we have 

lim / \Vv n (y)\ 2 dy = oo. 

Lemma 2.2. Let x E Q be an arbitrary point and a be a needle with tip at x. Let 
£ = {v n } be an arbitrary needle sequence for (x,a). Then for any point z G a and open 
ball B centered at z we have 

lim / \Vv n (y)\ 2 dy = oo. 

n^oo J Bnn 

From these lemmas we know that one can recover full knowledge of the given needle as the 
set of all points where the needle sequence for the needle blows up. This means that the 
needle is realized as a special sequence of harmonic functions without loosing information 
about the geometry of the needle. This is the new point added on the probe method in 
[18]. 

Definition 2.3. Given x G f2, needle a with tip x and needle sequence £ = {v n } for 
(x, a) define 

I(x,a,0n= [ {(A -A D )J n }f n dS,n=l,2,--- 
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where 

f n (y) =v n (y), y e dtt. 

{I(x, a, ^)n}n=i,2,- is a sequence depending on £ and a G JV X . We call the sequence the 
indicator sequence. 

The behaviour of the indicator sequence has two sides. One side is closely related to 
the function defined below. 

Definition 2.4. The indicator function I is defined by the formula 

I(x)= ( \VG(y-x)\ 2 dy+ I _\Vw x \ 2 dy, xett\D 
Jd Jq\d 

where w x is the unique weak solution of the problem: 

Aw = in Q \ D, 
dw d 

_ = -_«?(. -,) )<m I>, 
w = on<9f2. 

The function w x is called the reflected solution by D. 
The following theorem says that 

• one can calculate the value of the indicator function at an arbitrary point outside 
the cavity from A — A^; 

• the indicator function can not be continued across the boundary of the cavity as a 
bounded function in the whole domain. 

Theorem 2. 3. A We have 

• (A.l) given x G Q\D and needle a with tip at x if a n D = 0, then for any needle 
sequence £ = {v n } for (x,cr) the sequence {I(x, a, £) n } converges to the indicator function 
I(x) 

• (A. 2) for each e > 

sup I(x) < oo 
dist(x, D) > e 

• (A. 3) for any point a G dD 

lim I(x) = oo. 

X ><2 

Since mathematically Theorem 2. 3. A is enough for establishing a reconstruction for- 
mula of the cavities, in the previous applications of the probe method [10] we did not 
consider the the natural question 

• what happens on the indicator sequence when the tip of the needle is just located 
on the boundary of cavities, inside or passing through the cavities? 

However, in practice the tip of the needle can not move forward with infinitely small step 
and therefore in the scanning process with needle there is a possibility of skipping the 
unknown boundary of cavities, entering inside or passing through the cavities. So for 
the practical use of the probe method we have to clarify the behaviour of the indicator 
function in those cases. The answer to this question is 
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Theorem 2.3.B Let x eVL and a e N x . If x eVl\D and a n D ^ or x E D, then for 
any needle sequence £ = {v n } for (x,a) we have lim n >0O I(x, a, £)„ = oo. 

These two theorems are essentially a special case of the results in [18]. Finally we note 
that, as a corollary of Theorems 2.3 A and B one gets a characterization of D. 

Corollary 2.4 A point x G Q belongs to Q\D if and only if there exists a needle a with 
tip at x and needle sequence £ for (x,o~) such that the indicator sequence {I(x,o~,£) n } is 
bounded from above. 

Needless to say, this automatically gives a uniqueness theorem, too. 



3 Main idea. A new role of a Carleman function 

We think that the reader now understand the basic of the probe method and the role of 
the needle sequence. So the next problem is how to construct the needle sequence. 

There are several points for the meaning of the "construction" . It is well known that one 
can construct the needle sequence by solving infinitely many first kind integral equations 
in the sense of minimum norm solutions. The minimum norm solutions are given by a 
combination of the Tikhonov regularization method and Morozov discrepancy principle(see 
[25] for these concepts). However, from our point of view this does not yield an explicit 
needle sequence. One can also point out that the determination of the regularization 
parameter via the Morozov discrepancy principle itself is a nonlinear problem and seems 
impossible to find the parameter explicitly. 

So how can one construct the needle sequence explicitly? In this and following sections 
we always consider a geometrically simplest needle described in 

Definition 3.1 A needle with tip at x is called a straight needle with tip at x directed to 
uj if the needle is given by l x {uS) H Vt where 

I x {uj) = {x + tuj I < t < oo}. 

In two-dimensions in [22] we constructed an explicit needle sequence for G given by 

G(y) = (3.1) 
yi + m 

Note that this G satisfies (2.1) for any finite cone V with vertex at 0. 

For the construction we made use of a well known idea in complex function theory: 

1 1 



y — x (y — x') — (x — x') 
1 

, / x-x'\ 
(y-x)il- 



E 



y — x 
(x - x') m 

o (y - x') m+l 
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provided x' ~ x and x' is located on a needle with tip at x. Note that we used the same 
symbol for the complex number Z\ + iz 2 corresponding to a point z — (z±, z 2 ) T . 
However, in general we have to do this type procedure many times and therefore the 
resulted sequence is quite involved. 

In this section we give a different and extremely simple needle sequence and explain 
the idea behind. 

Definition 3.2. Given a unit vector uj = (oji, uj 2 ) t in R 2 and a e]0, 1] define 

E a (ryuJ)-l 

v{y; a, r, uj) = , r > 0. 

V 

This is a harmonic function in the whole plane. Note also that y ■ {00 + iuj ± ) = yuJ where 
UJ L = (-uj 2 ,uji). 

Theorem 3.1 Let x G Q and a be a straight needle with tip at x directed to uj. The 
sequence {v( ■ —x; a n , r n , uj) |q} is a needle sequence for (x, a) where a n and r n are suitably 
chosen sequences and satisfy 

• < a n < 1, r n > 

• a n — >■ and r n — > 00 as n — > 00. 

Proof. For each fixed a we have, as r — ► 00 

v( - — x; a, r, uj) — >■ G{ • — x) 



uniformly for y in any compact subset of R 2 \ C x (ira/2,uj). This is a direct consequence 
of the well known asymptotic behaviour of Mittag-Lefner's function [4, 8]: as \z\ — >■ 00 
with na/2 + e < |arg^| < n 



EJz) —7r(^aj + {w) (X2) 

where a and e are fixed and satisfies < a < 1 and < e < n — na/2. 

The choice of {a n } and {r n } depend on that of an exhaustion for Q\o~. More precisely 
let {O n } be a sequence of open subsets of Q\cr such that: for all O n C O n+ i, U n O n = Q\a. 
We call this sequence an exhaustion for Q \ a. Let {e n } be an arbitrary sequence of 
positive numbers such that e n — ► as n — > 00. The choice of {a n } and {r n } can 
be done as follows. Since 0\ C fl \ a, one can choose a small a e]0, 1[ in such a way 
that 0\ is contained in the set R 2 \ C x (na/2, uj). Set ai = a. Since v(y — x;ai,T,cu) 
converges to G(y — x) in H 1 (Q 1 ) as r — > 00, one can find a large r in such a way that 
\\v ( • — x; aii, t,uj) — G{ - — x) ||iji(Oi) < e i- Set T\ = r. For a 2 , • ■ ■ and r 2 , • • • we do a 
similar procedure. Then the obtained sequence satisfies the desired property. 
□ 

We see that v(y — x; a n , r n , uj) takes the special value at the tip x of the needle: 

v(y-x;a n ,T n ,uj)\ y=x = - TnU (3.3) 

r(i + a n ) 



and thus this yields 



v(y - x;a n ,T n ,uj)\ y=x _ 
hm — = — uj. 
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Therefore the leading term of the asymptotic behaviour of this special needle sequence at 
the tip x of needle a does not depend on the choice of {«„} and is uniquely determined 
by {r„}. 

However, the point of this section is not the statement of Theorem 3.1 itself. How 
could the author find this simple form? This is the main point. For simplicity let x — 0. 
First the function G(y) given by (3.1) is a fundamental solution for the operator d (ig- 
noring multiplying a constant, hereafter same). Since E a (0) = 1, the function 

E a (ryuJ) 



y 



(3.4) 



is also a fundamental solution for operator d. The point is: this fundamental solution 
becomes small outside a sector as r — > oo ([4, 8]). Since we can write 

Eajryu) = E a (TyuJ)) - 1 

y y y 

this means that the outside the sector 

E a (ryuj) - 1 



Note that this righthand side is just coming from the regular part of (3.4). 

The function (3.4) has been appeared in Yarmukhamedov's work [31] to give an explicit 
formula of the value of the solution of the Cauchy problem in a special domain for d. Such 
type fundamental solution for some operator is called the Carleman function (for the op- 
erator) and the formula is called a Carleman-type formula. See also [1] for Carleman-type 
formulae in complex function theory and [2, 15] for the stationary Schrodinger equation. 

So now the principle of the choice of needle sequence is clear. In particular, in three 
dimensions the principle is: choose (—1) times the regular part of a suitable Carleman 
function (for the Laplace operator). 

This principle is the conclusion of this paper and gives a new role for the Carleman 
function. We found a direct link between the Carleman function and the probe method. 
Needless to say, it is different from the idea in existing author's previous application [14] 
of the Carleman function to an inverse boundary value problem in unbounded domain. 
In the next section we consider the three dimensional case more precisely. 

4 Yarmukhamedov's fundamental solution and nee- 
dle sequence in three dimensions 

The following is taken from [30], however, therein the proof is not given. For reader's 
convenience we gave a direct proof of this theorem in Appendix. 

Theorem 4.1 Let A > 0. Let K(w) be an entire function such that 

• K(w) is real for real w 

• K(0) = 1 

• for each R > and m = 0, 1, 2 

sup \K {m \w)\e x \ Imw \ < oo. (4.1) 
\Rew\<R 



8 



Define 

2 <I>iri = , 

2 Jo W , /U/12 +u 2 



-2tt 2 $(x) = - r /m (^M) e "" + e XU du . (4.2) 



where w = x 3 + |x'| 2 + -u 2 and x' = (xi, x 2 ). Then one has the expression 

e iX\x\ 

$(x) = — — r + if (x) 
47r|x| 

where H is C 2 in the whole space and satisfies 

AH(x) + X 2 H(x) = mR 3 

and therefore $ satisfies 

A$(x) + A 2 $(x) + 5(x) = znR 3 . 

Hereafter we set $ = $x to denote the dependence on K and consider only the case 
when A = 0. In this case (3.2) ensures that K(w) = E a (rw) with r > satisfies (4.1). 
Yarmukhamedov [32] established that the function $x for this K with a fixed a is a 
Carleman function (for the Laplace equation). 

Definition 4.1 Given two unit vectors i?i and t? 2 in three dimensions and a e]0, 1] define 
v(y; a, r, tfi, 2 ) = - S^ K (y ■ &i, y &2, V (^i x t? 2 )) - - r > 

where K(w) = E a (rw). From Theorem 4.1 for A = one knows that this function of y is 
nothing but the regular part of $x and harmonic in the whole space. 
In this section we prove 

Theorem 4.2 Let x e VL and a be a straight needle with tip at x directed to uj = i?x x 
Then the sequence {v( ■ — x; a n , r n , i? 2 )|n} is a needle sequence for (x, a) with 

G(y) 



4:7c\y\ 



where a n and r n are suitably chosen sequences and satisfy 

• < a n < 1, r n > 

• a n — >■ and r n — >■ oo. 

Proof. 

First we prove that, for each fixed a as r — >■ oo 

v( - — x; a, t, i? 2 ) — >■ G( • — x) 



uniformly for y in any compact subset of R 3 \ C x (na/2, uj). 

Given r\ e]0, 7r[, r > let 7(77, r) denote the contour that originates at 00, runs 
in toward the origin just above the half line z = pe~ tr, ,p > r, the part of the circle 
z = re ld , \9\ < i] the origin counterclockwise, and then returns to 00 just above the 
half line z = pe tr, ,p > r. The contour 7 splits the complex plane into the two simply 
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connected infinite domains D~(r),r) and D + (i],r) lying respectively to the left and the 
right of 7(77, r). 

From Mittag-Lefner's integral representation (p. 206 of [4]) and analytic continuation, 
it follows that, for any 77 e]7r/2, ir[ 

Let r > 0. We see that: if z e D~(ar],r a ) \ B ra (0), then tz e D~(ar],r a ); if z G 
D~(arj, r a ), then so is z. And it should be pointed out that E a (z) is real for real 2. From 
these we have, for all z E D~(arj,r a ) \ _B r «(0), 

1 (E a (rz)\ _ 1 /■ (C Q -2rRe2)C a " 1 e^ 



Imz \ Z J 2ixi\z\ 2 Jj(r,,r) (( a - Tz){C, a - TZ) 

Here we restrict the interval where r\ belongs to: r\ <E]n/2, ir/(2a)[. Let 6 satisfy 
ar] <6 < tt/2. Then we have, for all w G D~(9, r a ) \ B r <*(0) and r > 1 

inf \C-tw\>Ct\w\ (4.5) 

C€7(»7,r) 

where 

C = min {1 , sin {9 — ar])}. 

T 

This can be proved as follows. It suffices to consider the case when 9 < axgw = tp < ir. 
We divide the case into two subcases. First let 9 < ip < ai] + tt/2. In this case the circle 
\z — w\ = \w\ sm(ip — ar)) has a single common point with the half line z = ke tav , r a < 
k < 00. Thus this yields 

inf \( a — w\ = \w\ sin (ip — ar]). 

Since 9 — ai] < ip — at] < tt/2, we obtain 

inf \( a — w\ > \w\ sin (9 — ar]). (4.6) 

Next consider the case ai] + n/2 < ip < ir. A simple geometrical observation yields 

inf \( a -w\ = \w-r a e iav \>\w\-r a . (4.7) 

Now let t > 1. Then rw still belongs to D~(9,r a ) \ B r a(0) and having the same arg as 
w. Therefore inserting rw instead of w into (4.7) we obtain 

inf \C -tw\ > t\w\ -r a > (r- l)\w\. (4.8) 

C€7(»?,r) 

Now a combination of (4.6) and (4.8) gives (4.5). 
Now let y satisfy 

(y — x) • u < \y — x\ cos 9 (4.9) 

and 

r a < \y - x\ < R a . (4.10) 
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Set w = (y — x) -a; -My \(y — x) • -#i| 2 + \ (y — x) ■ -# 2 | 2 + u 2 . This w belongs to D (6,r a ) \ 
B r c(0). Therefore from (4.4) and (4.5) we obtain 

1 1 i m fg#^i< J; , / (icr+2i2-)ici- i e Ite C| dC |. 



Imw \ w J 2ir r (r 2a + u 2 ) 2 Jjfar) 

The integral in this right hand side is absolutely convergent since n/2 < r/ < n. Therefore 
we have , 

$ K ((y - x) ■ (y-x)- # 2 , (y-x)-u) = O(-) 

r 

as t — > oo and the estimate is uniform for y — x with (4.9) and (4.10) for each fixed 6>, 
f], r, R and a. 

It is easy to see that, choosing suitable 6, rj, r, a large fixed R and a, one can generate 
an exhaustion of the domain Q \ a. Therefore from the principle explained in Section 3 
we conclude that the sequence {v ( • — x; a n , r n , ^2) |o} converges to G(y) in L 2 Qc (f2 \ 
a) for suitably chosen sequences {a n } and {r n }. Then a standard argument gives the 
convergence in H\ oc (Q \ 
□ 

Remark 4.1. The proof presented here can be considered as a minor modification of that 
of Lemma in [32]. Anyway the new point that should be emphasized is the relationship 
between the construction of a needle sequence and a Carleman function. 
Definition 4.2. We call the needle sequence {v ( • — x; a n , r n , $1, "$2)!^} given in Theorem 
4.2 a standard needle sequence for (x, a) for the Laplace equation. 
Choosing K{w) = 1 in § K , we have 

*i(v) = T~\ 1 • 
47r|y| 

This gives the expression: 

~ 9 / n /*°° T (EJtw)-1\ du , t „„. 

2* 2 v(y;a,T,# 1 ,& 2 = / Im ^ - 4.11 

7o V w / v /|y^i| 2 + |ytf 2 | 2 + w 2 

where 



w = y • u + i J\y ■ $i| 2 + \y ■ $ 2 \ 2 + u 2 . 



The advantage of the standard needle sequence is that one can exactly know the precise 
values in a closed-form on the needle and its linear extension. 

Theorem 4.3 Let u = #1 x $ 2 . We have: 
(1) if y — x + sou with s ^ 0, then 

v(y - x- a, r, ^, 2 ) = E ^ ~ 1 ■ ( 4 .12) 

47T S 



(%) if y = x, then 



v(y-x; a, r, ti x , fl 2 ) | y=3! = ^ - . (4.13) 
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Proof. 

From (4.11) we have the expression of v(y — x; a, r, $1, i? 2 ) i n both cases (1) and (2): 

« 2 / n nx r, f E a (T(x 3 +IU)) - 1\ rf« 

27T 2 w (y-x;a,r,i?i,i?2 = / Im — — (4.14) 

Jo \ x 3 + iu J u 

here we set x 3 = (y — x) • uj for simplicity of description. One can write 

/■°° Im / ff Q (r(:r 3 + m))-l \ da = l_ ^ / /— + P/A E a (r(rr 3 + Q) ~ 1 (4 15) 

First consider the case when x 3 ^ 0. Then ( = —x 3 is a removable singular point for the 
integrand. Therefore, from the Cauchy's integral theorem and (3.2) we see that the right 
hand side of (4.15) becomes 

1 Um f E a (r(x 3 + £))-! n E a (rx 3 ) - 1 



2l e|0 j£=ee ie , |0-7r|<7r/2 (x 3 + ()( " 2 :T 3 

This together with (4.14) and (y — x) ■ uj = s for y = x + suj yields (4.12). 

A similar computation in the case when x 3 = yields (4.13). 

□ 

The formula (4.13) corresponds to (3.3). 

Remark 4.2. From (4.12), (4.13) and the trivial estimate 

x n 

E a {x) - 1 > — — -, x > 0, n = 1, 2, • • • 

1 (1 + an) 

we see that the standard needle sequence on the needle blows up in a point wise sense: 

for all points y on l x (uS) (see Defintion 3.1) lim„ >oc v(y — x; a n , r n , $2) = +00. 

Remark 4.3. We can know also the values of Vv(y — x : a, ti?i, #2) on the line y = 
i + sw (—00 < s < 00) with a; = i?i x i? 2 . Define 



_ , s /•°° T / ^^(s + iVp + w 2 )) \ *x 

h(p,s)= Im V \ . p > -1, -00 < s < 00 



and 



J 2 p, s = / Im , • r^ r~^ r^i > < p < 1, -00 < s < 00. 

jo \ s + 2 VP + M J VP + U 

Using (3.2) and the asymptotic expansions of the derivatives of Mittag-Lefher's function 
outside the sector |argz| < na/2, we conclude that h(p, s) is a smooth function of (p, s). 
Following the argument in Appendix A. 2 we have 



>(P> S ) = - / 

JO 



du 



P + S + M 

-n+1 



^r(l + a(n + 1)) ^ odd j\(n-j)\ 1 ; io ^ J 
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Since the second term is smooth for (p, s) G R 2 , this yields that the function 

du 



I 3 (p,s) = I 2 (p,s) + I 
Jo 



p + s 2 + u 2 



has a smooth extension to R 2 . Since 

1 \ du r°° du 



r Im , \"=-r 

Jo \ s + + u z I v p + u z Jo 



S + iy/p + U 2 ) \/p+U 2 JO p + S 2 +U 2 ' 

we have the decomposition of (4.11): 

2n 2 v(y, «, #i, # 2 ) = h(p, s) + I 3 (p, s) (4.16) 
where p = \y ■ -#i| 2 + \y ■ $2| 2 and s — y • uj. Therefore we obtain 

2ix 2 Vv{ y] a, r, # u # 2 ) = —(A + J 3 ) + 2— (A + J 3 ) {(y • W + (y ■ d 2 )§ 2 }. 

This together with (4.12) and (4.16) yields that for y = x + su with s ^ 

Y7 ( 9 ^ d j E a (rs)-l \ 

Vv(y-x;a,T,AM = Ts { ^ jo;. 

Moreover, letting s — ► 0, we have 



r 2 



V%-x;a,^ 1 ,^)|^ = M(i + 2a) ,. 

Therefore Vt> on the line y = x + s u (— oo < s < oo) is parallel to the direction of the 
line. In particular, using the power series expansion of Mittag-Leffler's function, we see 

that Vv{y — x; a, r, i? 2 ) • w > and lim n >oc Vi>(y — x; a n , r n , #2) • ^ — +°° on the 

set I x {uj). 



5 Needle sequence for Helmholtz equation and Vekua 
transform 

The needle sequence for the Helmoltz equation can be defined as same as the Laplace 
equation. 

Let G(y) be a solution of the Helmholtz equation Av + X 2 v = in R 3 \ {0} such that, for 
any finite cone V with vertex at G satisfies (2.1). Hereafter we fix G. 
Definition 5.1. Let er G N x . We call the sequence £ = {v n } of if 1 (fi) solutions of the 
Helmholtz equation Av + \ 2 v = a needle sequence for (x, a) for the Helmholtz equation 
if it satisfies, for any compact set K of R 3 with K C \ a 

JimJIM -)-G(--x) \\ LHK) + ||VR( • ) - G( ■ - x)}\\ LHK) ) = 0. 

Remark 5.1. Any needle sequences for the Helmholtz equation with an arbitrary A blows 
up on the needle in the sense of Lemmas 2.1 and 2.2. See [18] for the proof. 
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The existence of the needle sequence for the Helmholtz equation has been ensured under 
the additional condition on A (see Appendixes of [11, 18] for the proof): 

• A 2 is not a Dirichlet eigenvalue of —A in Q. 

However, the proof therein does not give us any explicit form of a needle sequence. Since 
the needle sequence for the Helmholtz equation plays a central role in the probe method 
applied to inverse obstacle scattering problems (see [11, 18] again and [21] for a recent 
application), it is quite important to obtain an explicit one. 

In this section we give explicit needle sequences for the Helmholtz equation for all 
straight needles and all A(> 0). There is no additional condition on A. For the construction 
we do not make use of Theorem 4.1. Mittag-LefHer's function with a < 1 does not satisfy 
the condition (4.1) when A > 0. Therefore one can not substitute K(w) = E a (rw) 
into (4.2) to construct a needle sequence for the Helmoltz equation Au + X 2 u = in 
three dimensions. In this section, instead of (4.2) we employ an idea of making use of a 
transformation introduced by Vekua [26, 27]. 

Let A > 0. The Vekua transform v i — > T\v in three dimensions takes the form 

Txv(y) = v(y) - ^ f^v^J^X^y/T^t) ^^dt 

where J\ stands for the Bessel function of order 1. The important property of this trans- 
form is: if v is harmonic in the whole space, then T\v is a solution of the Helmholtz 
equation Au + X 2 u = in the whole space. 

Definition 5.2. Given two unit vectors i?i and $2 in three dimensions and a G]0, 1] 
define 

v\y\ a, r, A, #2) = T x v(y; a, r, & u tf 2 ), r > 0. 

This function of y satisfies the Helmholtz equation Am + X 2 u = in the whole space. 

Theorem 5.1 Let x G Q and a be a straight needle with tip at x directed to uj = $1 x i? 2 - 
Then the sequence {v x ( - — x; a n , T n , $1, $2)!^} is a needle sequence for (x,a) for the 
Helmholtz equation with G = G\ given by 

where a n and r n are suitably chosen sequences and satisfy 

• < a n < I, T n > 

• a n — ► and r n — > 00. 

Proof. 

It suffices to prove the theorem in the case when x — 0, $1 — (1,0, 0) T , $2 = (0, 1,0) T 
and thus u = (0, 0, 1) T . 

First we show that for each fixed a as r — > 00 the function v x (y; a, r, i? 2 ) converges 
to the function v(y) defined by 

v(y) = ^(y) - ^ jf 1 9 1 ( t y)J 1 (X\y\y/T^T)^L.dt, y ? (5.1) 
uniformly for y in any compact subset of R 3 \ Co(na/2, uj). 
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From Theorem 4.2 we have, <&k{u) with K(z) = E a (r z) converges to uniformly for y 
in any compact subset of R 3 \ Co(na/2, ui) as r — > oo. Since $i(y) — v(y; a, r, $2) = 
$k(u), it suffices to prove that the following function of y 



R(y;a,r) = j\ K (ty) J^\\ y \y/T^t) J^dt, y ± 0, (5.2) 



converges to uniformly for y in any compact subset of R 3 \ Co(vro;/2, a;). 

Let 77 G]7r/2,7r[. Let r and i? be arbitrary positive numbers with r < R. Let 9 G 
]a?7, 7r/2[. Let y satisfy (4.9) and (4.10) with x = and uo = (0, 0, 1) T . Then the non zero 
complex number 

w = y 3 + i\jyl + yi + u 2 , u > 

satisfies 9 < |argu>| < ir. Then we see that, for all r > tw G D(ar/,r a ). Thus (4.3) 
gives 

EJ TW ) = — 4 — d C- 



Since E a (z) = E a (z) for all z, we obtain 



jo 1m w \ w J 
' 1 I (C;-2rRe,0C°-' e < 1 

•Avfa.r) £ Q -TIC -Tl 



27ri7o jJH 2 AOvO (C° -™)(C a -rw) 
Let < i < 1 . A change of variable u — > tu yields 

-2. = -_ y o {-gr^j^ {(a _ Ttw){(a _ Ttw) dC] du. (5.3) 

Divide 7(77, r) as 71 + 72 + 73 where 

71 : C — se " ') t < s < 00, 
—72 : C — se~ ir? , r < s < 00, 

73 : C = re^, |<p| < 77. 
Then from (5.2) and (5.3) we have the expression 

4ir 3 iR(y; a n , r„) = i?i + i? 2 + #3 

where 

R i= r 1 (r{' 3 r /"-^f'^^/'^f'^j. (5.4) 

J io [Jo \\y\ 2 + u 2 J^(C a -rtw)(C a -rtw) S J ^/^T^) / 

From a simple geometrical observation, for all t G [0, 1] and u > we have the estimates 
|e ia?J — tw\ > sin (0 — arj) and |e im? — tw\ > sin (0 — arj). Note that the common right hand 
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side of theses estimates gives the distance between the point e iari and the set {w | 9 < 
\aigw\ < 7r}. 

Then these estimates together with a change of variable yield 

(C Q - 2Tty 3 )( a - 1 e<; f°° s a - 1 e scosr >( s a + 2T\y 3 \)ds 



f {Q"-^Tty 3 )Q—e" < r 

Ai (C a -Ttw)(C a -TtW) ~ Jr 

1 r a 

a Jr a 



hi (C° - Ttw)(( a - rtw) Jr \s a e iar > -Ttw\\s a e iar > -rtw\ 

1 roo e Tl/ ^ 1/acos ^(e + 2|y 3 |)^ 



It \e tar i - tw\\e iar i - tw\ 
a sin (9 — arj) Jo 

I roo 

< T7~ n \ — Ti~ / e scos "(r- 1/Q s + 2R a )ds = 0{T- l ' a ). 

~ a sin 2 (6» - ar])r 1 / a Jo 

Similarly we have the same estimate for the integral of the same function on 72. Therefore, 
using the well known estimate 

\Ji(s)\<^,s>0, 

we obtain, for j — 1,2 Rj = 0(t _1 / q ) as r — > 00. This is a uniform estimate with y 
satisfying (4.9) and (4.10) for x = and cu = (0, 0, 1) T . 

Next we show that R 3 = 0(r -1 / 2 ) as r — > 00. 
We have 

(C-2r^ 3 )C- 1 e f 



I(w,Tt) = \ -dCI 



<r a e 

J-T) 



;r cos <(j 



jrV^ - 2rty 3 \dip 



dip 



e ia<p _ T f w \\ r a e ia V - T fyj\ 



<(rRni + 2rt)£ r — 
and thus this yields 

/-i / r°° I{w,rt) du \ IMXlylVT^t^dt 

/ n x«Ai? a /" 1 l + 2rt , /■«> du ( rv dp 
< (rR) a / i^dtl -- r / -. : rr — = r 

~ 2 Jo v 7 ^ -A) y 2 + m 2 \J-v \r a e ia ^ - r i w \\ r a e ia <f - r tw 



(5.5) 
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A change of variable gives 



1 1 + 2r t , f°° du ( rv dip 
at 



y/i Jo \y\ 2 + u 2 \J- V \r a e ia *> - t t w\\r a e ia *> -rtw\ 



1 n 
\/t Jo 



1 + 2£ f°° du ( fv dtp 
-aq 



Vt Jo v 7 ^ Jo \y\ 2 + u 2 \J-v \r a e ia v - q~w\\r a e iai ? - q~w\) 

< 1 [°° 1 + 2 ^ d£ [°° dU ( r ^ \ 1 (I I 

~ Vr Jo v 7 ^ Jo \y\ 2 + u 2 \J-n \r a e ia ^ - ^w\\r a e ia ^ - ) ~ ^/r 1+ 2 ' 



where 

1 + 2£ f°° du ( n dip 



(5.6) 



11 J2 ^ Jo \y\ 2 + u 2 \J- V \r a e iai ? - q~w\\r a e iai ? - £w\ 

and 

I r 2 1 + 2^ r°° du ( rv dip 

2 Jo \^ Jo \y\ 2 + u 2 \J-r, \r a e iai ? - 



£ w\ \r a e %olv — £ w\ 



If £ > 2, then we have 



\r a e iaip -£w\> £\w\ -r a > £\w\(l - -jA) > £\w\(l - -A) > — . 
Therefore we get 

roo \ _|_ 2£ r°° du 

On the other hand, a geometrical observation gives, for all £ > 

\r a e iav -H > r a sin(^-«r7). 

This right hand side is just the distance between two sets {w \ 9 < |argiu| < n} and 73 in 
the complex plane. Therefore we obtain 

2rj f 2 1 + 2£ f°° du 

h ~ r 2 « sin 2 (9 - or,) Jo ~7T ^ Jo ^+u~ 2 < °°' ^ 

Now from (5.5) to (5.8) we obtain the desired estimate for R 3 . 

Summing up, we have R = OiT' 1 ' 2 ) as r — > 00 and this is a uniform estimate 
with y satisfying (4.9) and (4.10) for x = and u = (0,0, 1) T . Therefore the function 
v A ( • ; a, r, $2) converges to the function v(y) uniformly for y in any compact subset 
of R 3 \ Co(na/2, uj). Since a can be arbitrary small and the function v x ( ■ ; a, r, #1, i? 2 ) 
satisfies the elliptic equation Am + A 2 w = 0, a standard argument yields that, for suitably 
chosen sequences {«„} and {r„} the sequence v A ( • ; a n , r n , i? 2 ) converges to u(y) in 

The next point of the proof is that the right hand side of (5.1) coincides with G\(y). 
This is a direct computation. The second term of the right hand side of (5.1) becomes 

A f 1 MMy - x\y/T=1) dt = A f\l - w 2 )' 1 ' 2 J^y - x\w)dw. (5.9) 
8tt Jo _ t ) 4tt Jo 
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Here we know that 



s [\l - w 2 )- 1/2 J 1 (sw)dw = 1 - cos s. (5.10) 
Jo 

This is an implication of (30) for (5 = 0, a = s and v — 1 on p. 337 of [3]: 

jy-xr i/2 J^)dx = \{j 1/2 fyY 

where Ji/ 2 stands for the Bessel function of order 1/2 and has the explicit form 



Jiniw) = \ sin w, w > 0. 

V ttw 

Using (5.1), (5.9) and (5.10), we conclude that v (y) = G\(y). 
□ 

Definition 5.3. We call the needle sequence {v A ( • —x; a n , r n , $2)!^} given in Theorem 
5.1 a standard needle sequence for (x, a) for the Helmholtz equation Au + \ 2 u = 0. 

Therefore now we can say that the probe method applied to the inverse obstacle 
scattering problems becomes a completely explicit method if one uses only straight needles 
and standard needle sequences. In this case this explicit method gives information of the 
location and shape of unknown obstacles more than the convex hull. 

Some remarks are in order. 
Remark 5.2. From Theorem 4.3 and (5.10) we get the precise values of v x (y—x; a, r, #1, i? 2 ) 
on the line y = x + su(—oo < s < 00) with uj = $1 x $ 2 - The results are: 

if y = x + s uj with s 7^ 0, then 

x ( a a x 1 E a (rs) - cos As 
v {y-x; a, r, tfi, $ 2 = 

47T S 

f\l - w 2 )- 1/2 E a (T(l - w 2 )s)J 1 (Xsw)dw; 
4tt Jo 



if y = x, then 



v x (y - x;a,r,i?i,T? 2 ) 



r 



|J/=X 



47rr(l + a) 



Moreover, from (4.16) and the fact that Ji(s) is an odd function we see that Vv x (y — 
x; a, r, i? 2 ) 011 the line y = x + s (—00 < s < 00) is also parallel to uj. In particular, 
we have 



Vv (y -x; a,T,$i,$ 2 )\y=x 



T 



2 



4vrr(l + 2a) 
Remark 5.3. Since the function 



LO. 



I 



47r|y| 

satisfies the Helmholtz equation in the whole space, we see that the sequence {v x (y — 
x\ a n , T n , i? 2 )} given by 

w (y - x;a n ,T n ,^i,^ 2 ) +«— — -, n=l,---,ye$2 

47r \y — x 
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is also a needle sequence for (x, a) with G — G\ given by 



Gx(y) = 



47r|y| 



This G\ is nothing but the fundamental solution for the Helmholtz equation and sat- 
isfies (2.1) for any finite cone V with vertex at 0. We call also the sequence {v x (y — 
x; a n , T n , ^2)} the standard needle sequence for (x, a). 

Remark 5.4. Fix a G ]0, 1[. As a corollary of the proof of Theorem 5.1 we see that: the 
function 

*{y)=G x {y)-v\y;a,T,-& u -& 2 ), y^O 

with a large parameter r is a Carleman function (for the Helmholtz equation) for a special 
domain whose boundary consists of a part of the conic surface and a smooth surface lying 
inside the cone of axis direction i?! x i? 2 an d aperture angle 7ra. This is an extension of a 
result in [32] to the Helmholtz equation and new. As pointed out in the beginning of this 
section Theorem 4.1 for A > does not give this result. 

Remark 5.5. For a recent application of the Vekua transform to a boundary value 
problem for the Helmholtz equation see [6]. In [19] we applied the idea of the Vekua 
transform to the construction of the Herglotz wave function that approximates in a disc a 
complex exponential solution of the Helmholtz equation. This is an application to inverse 
obstacle scattering problems in two dimensions. 



6 Conclusion and open problems 

We could find an explicit needle sequence for a needle given by a segment in three dimen- 
sions. This makes the probe method completely explicit as same as the enclosure method. 
It is based on a principle that connects the regular part of a Carleman function with a 
needle sequence in the probe method. 

Some further problems to be solved are in order. 

• A mathematically interesting question is: can one give an explicit needle sequence having 
a closed-form for a general needle not necessary given by a segment in three dimensions? 

• Using the standard needle sequence for the Helmholtz equation Au + X 2 u = 0, one can 
construct the corresponding indicator sequence for some inverse boundary value problems 
that are the reduction of inverse obstacle scattering problems at a fixed wave number 
to a bounded domain. Can one prove the blowing up of the indicator sequence without 
assuming that A is small when the tip of the needle is located on the boundary, inside or 
passing through unknown obstacles. For the case when A is small see [18, 21]. 

• Construct a needle sequence for the stationary Schrodinger equation — Au + V(x)u = 
or the equation V • 7(2;) Vw = as explicit as possible. This is important for the recon- 
struction problem of discontinuity embedded in an inhomogeneous conductive medium. 
See [16] for the probe method applied to this problem. 

• Do the numerical implementation of the probe method by using the standard needle 
sequences and compare the results with the previous numerical study done in [7, 9] of the 
probe method in two dimensions. 
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A Appendix. Proof of Theorem 4.1 

A.l Preliminaries 

Given an entire function K(w) set 

'> = / 

Jo 

where \& is a smooth function on [0, oo[ and 



, K(w) tt(Au) , 
®(x) = I — — . — — du 

w -nx'\ 2 + u 2 



w = X3 + iy \x'\ 2 + u 2 . 

Set 

f(w) = . 

w 

Here we show that, under suitable growth conditions for f(w) as u — > 00 and the choice 
of fy(\u) the function satisfies 

A$ + A 2 $ = 

for x' ^ 0. 

First we do the following formal computations. We have 

**=r/» nxu) du. (A.i) 



dx\ jo /i r /,v + // 



Let k = 1,2. We have 

|i= f f( w )i Xfc - *(\u)du- I" f(w) . Xk *(\u)du. 
dx k Jo v 1 \x'\ 2 + u 2 v 1 Jo v (J\x'\ 2 + u 2 f 

Then we obtain 

<9 2 $ t°° f " t x l rt\ \ a f°° f>( -s-^(^ u )du 



2 + u 2 



0'® roc x * roc w ( 

*l = -i /(tt,) ^FT^* (Au) *' + i /W 'R 

/•oo 2t^ roc <y»2 

/ /12 ■ 2 x 2 ^(Ati)dti- / f»i 1 ^ (Xu)du 

Jo If I + u A y Jo If I + u A y 

f°° , V(\u)du f°° ,. , xl T/ . 
/ /(w) ^_ ' + / /( w )3— == fc m(\u)du. 

Jo (^\x>\ 2 + u 2 f Jo (^\x'\ 2 + u 2 f 
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This yields 

2 q2$ f°° s \x'\ 2 r°° s V(\u)du 

f°° \x'\ 2 f°° \x'\ 2 

-2 f/H ^ Xu)du +3 f/H - W? *(\u)du 

Jo ^\ x >\2 + M 2)3 7 (^(x'P + M 2 ) 5 

/■oo 3 It' 1 2 9 

= / /( W ){ ; ' ' }tt(Au)rftt 

7o r./l7T'|2-Hi.2^5 r./lT'|2 -I- 71.2^3 



x'| 2 +W 2 ) 5 (y|x'| 2 + W 2 ) ; 

2 3 It' 1 2 

Z" 00 // la;'! 2 

- / / ( w ) ' 1 ^V(\u)du 

Jo ( y /\x'\ 2 + U 2 ) 3 

\ |x'| 2 -2M 2 ■ f°° h, ^ ^2u2 -\ x '\ 2 Tf , \ i 

f-OC It' I 2 

- / f (w) . 1 l ^^V(\u)du. 
Jo f. MIU^3 



|X'| 2 +M 2 ) 

From (A.l) and (A. 2) we have 



A$= / f( w ) ' 1 tt(A«)(fa + z / fH ., „. ' ' - ^{\u)du 







2 -2u 2 




1 2 +M 2 ) 5 



2w 2 - 


x'\ 


2 


(\x'\ 


2 + u 2 ) 2 



r°° ,, I t 

+ / / (w){ - 1 1 }9(\u)du 

Jo ^\x'\ 2 + u 2 {^\x'\ 2 + u 2 f 



r°° a ii 2 
+ / f (w) , ^^^(Xu)du. 

Jo (J\ X ,\2 + U 2Y 



(A.2) 



(A3) 

f°° r/ \ W\ 2 -^u 2 . , .2u 2 -\x'\ 2 
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Integration by parts yields 

rOC r.,2 rOO f) ni 

Jo v U^ + u 2 )* v J Jo du XJK Jf \x'\ 2 + u 2 v ; 

- ' r + a r ^ 

Note that we assumed that 

u 

lim /'(«))— »*(Au) = 0. 

In what follows we always assume this type growth condition for f(w) to delete a term 
coming from integration by parts. 
Therefore (A. 3) becomes 

roc I nJ 1 2 91/2 

A$= / /(w) 7 1 =^tt(A«)<fa 

JO r ./It'|2 4- 7,2^5 



x'| 2 + 1i 2 ) £ 



Integration by parts yields 



Vo /H (M 2 + , 2 ) 2 = io ^ {/H} (Vl^| 2 + . 2 )3 
r/lx'|2 _|_ 7/2^|3 Jo <9« ( /U/|2 , ,,2\3 



- f(w) r— ty(\u)du-\ / /(m) — v ' 

Jo (yj\x'\ 2 + v?) 5 J ^\ x ,\2 + u 2y 



du. 



Therefore we obtain 



A$ = -A/ /H ; du + iX f(w)-—^ — ^dw. 



x'\ 2 + U 2 )' 



2 + u 2 



Here we assume that 

tf'(0) = 0. (AA) 
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Then integration by parts yields 



= /H-7=^= g°- / /H — 

./V 2 + u 2 J o du ,/\x^ + u 



x'\ 2 + V? 
' ' =du 



J- u*'(AiQ du _ x roo tf, Xu) ^ 

Jo (yVp + w 2 ) 3 Jo ^\x'\ 2 + u 2 



and thus we obtain 

roo 

A$ = -A 2 / /H- 

= (A6) 
A$ + A 2 $ = 



A$ = -A 2 /°° /H 1 j =dtx. (A5) 
Therefore if \I/ satisfies (A. 4) and 
then from (A. 5) we obtain 



for x' 7^ 0. 

The ^ that satisfies (A. 4) and (A. 6) is given by 

C{e u + e~ u ) 

and thus we choose 

*(Au) = C(e Au + e- Au ) 

where C is a constant. Then it is easy to check that the formal computations given above 
can be justified under the growth condition (4.1). 

Remark A.l. We give some historical remarks. In [28] he considered the case when 
A = 0. Therein he gave a proof for the fact that $ is harmonic for x' 7^ 0. However, 
the proof uses a change of variables and can not cover the case when K(w) satisfies (4.1) 
with A = 0, for example, K(w) = e mw with m > as pointed out in Appendix A of [14]. 
However, in [29], he had given a direct proof in that case. Using the idea of the proof, 
one can easily obtain that, in the case when A = $ is harmonic for x' 7^ provided 

Vi?>0 sup \K (m) (w)\<oo 
\Rew\<R 

for m = 0,1,2. 

Except for the explanation for the choice of a suitable \& the proof presented in this 
subsection follows that of [29]. 

A. 2 Extracting singularity at x = 

We assume that K(w) is real for real w. 
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Set 



^/M 2 + w 2 

Here we extract the singularity of this function at x = 0. 
Since 

K(a) = £ ^ ro » 



n=0 



we have 



if(w) 



if(0) 



Z^o Til v 



X3 + iJ\x'\ 2 + U 2 



n=2 



K(0)(x 3 -iJ\x'\ 2 + u 2 ) 



\x\ 2 + u 2 
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All coefficients K^(0),n — 0, 1, • • • are real and this yields 
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Therefore we obtain 
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= -*(°)/r^*< 
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modulo a C°° function on the whole space. Since 
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Jo 
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is C 2 on the whole space provided if satisfies a suitable growth condition, we obtain 
that 
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du 
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modulo a C 2 function on the whole space. 
Now we take 

tf(Au) = (e Xu + e- Xu )/2. 

Let < \x\ < 1. Then 

tf(Au) , 1 r 1 e Xu , 



Jo X 2 + U 2 2 7-1 



d 2 + M 2 



= ™ ReS (^T^ )Ul *' ~ 2 Jz=e»,0<0<n JxYT^^ 
JX\x\ ■ r7T „A(cos e+i sin 6) „i\\x\ 

' / -<'"<W 



2 Id 2 io 



|rr| 2 + 2 |x| 

modulo a C°° function on |x| < 1. Therefore from (A. 7) one concludes that 

K(0)tt e iX W 



\x\ 



modulo a C 2 function on the whole space provided K(w) satisfies (4.1). In particular, if 
K(0) = 1, then 

1 _ g^W 

modulo a C 2 function on the whole domain. 

Using the results in subsections, one obtains Theorem 4.1. 
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